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Abstract. We investigate spin currents and spin-current correlations for double 
quantum dots in the spin blockade regime. By analysing the time evolution 
of the density matrix, we obtain the spin resolved currents and derive from a 
generating function an expression for the fluctuations and correlations. Both the 
charge current and the spin current turn out to be generally super-Poissonian. 
Moreover, we study the influence of ac fields acting upon the transported electrons. 
In particular, we focus on fields that cause spin rotation or photon-assisted 
tunnelling. 



1. Introduction 

The recently achieved access to individual electron states of quantum dots [1] spurred 
the interest in the control of single electron degrees of freedom in nanoconductors 
providing a designable alternative to atoms in quantum optics. In particular by 
applying ac fields, interesting effects such as coherent destruction of electron tunnelling 
by electric ac fields [2-7] or electron spin rotations by crossed magnetic ac and dc 
fields [8-10] have been reported. A double quantum dot with up to two electrons 
provides a perfect framework for spin manipulation via spin blockade [11-14]. These 
systems are designed such that the only state located in the transport window is the 
doubly occupied singlet state. In this particular configuration, for electrons occupying 
inter-dot triplet states, i.e. two-electron states in which the electrons have the same 
spin orientation, inter-dot tunnelling is suppressed due to the Pauli exclusion principle. 
Tunnelling current in ac electric fields [4,15-17] can be tuned by means of the frequency 
and intensity of the field, which yields interesting features like for instance charge 
localisation within the quantum dot structure (dynamical localisation) and suspension 
of Coulomb and spin blockade. 

Spin blockade also plays an important role for coherent spin manipulation, as 
demonstrated by electron spin resonance experiments where a resonant ac magnetic 
field — perpendicular to a dc one — acts on the double quantum dot [9] giving rise 
to coherent single-spin rotations. Under certain configurations, the singlet states 
are populated such that spin blockade is lifted and, thus, an electronic current 
flows through the system. Then the current exhibits signatures of an interplay 
between coherent spin rotations due to the ac magnetic field and coherent electron 
delocalisation due to interdot tunneUing [10]. 
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Figure 1. Double quantum dot in spin blockade configuration. The chemical 
potentials in the contacts are adjusted such that the doubly occupied intra-dot 
singlet state is the only one located in the transport window. As electrons 
fall in the inter-dot triplet subspace, inter-dot tunnelling is forbidden by spin 
conservation. Consequently, transport is suppressed. 



In this paper, we analyse the spin current and the spin current fluctuations, 
as weU as their correlations in double quantum dots in the spin blockade regime in 
the presence of ac electric and magnetic fields. We demonstrate that spin current 
correlations provide information about the spin dynamics. Moreover, we also consider 
the case of spin relaxation which stems from spin flip processes accompanied by the 
coupling to a bosonic heat bath. This heat bath together with the fermionic electron 
reservoirs forms an unconventional environment for the quantum dots. 

2. Model and method 

We consider two weakly coupled quantum dots connected to two fermionic leads (cf. 
figure [1]) described by the Hamiltonian 

H = Hq + iJLR + Ht + -ff leads, (1) 

where Hq — J2icr ^t'^'^la'^io- + Z^i t^j^»T"^ii + C^lr"^l«r describes the two isolated 
quantum dots, Hlr = ~ X]cr(^LR^Lcr^R'^ + H.c.) is the inter-dot coupling. Ht = 
Sje{L R}ktT(7i'^Lo-^'<T + H.c.) represents the tunnelling between the double quantum 
dot and the leads described by -fficads = J2ik<j '^ik'jdl^.^dika- £ia is the energy of an 
electron located in dot i with spin a, and Ui (C^lr) is the intra-dot (inter-dot) Coulomb 
repulsion. In order to allow for up to two electrons in the system (one in each dot), the 
chemical potentials of the leads, fii, must satisfy the condition et < C/lr < Si + Ui 
and Hi < Ei + 2Ulr- In addition, a dc magnetic field is introduced along the z- 
axis in order to break the spin degeneracy by a Zeeman splitting = giBz.i, i-e. 

Hb = J2i ^i^z^ where Si = i X]<t<t' '^lcr'^<^<^"^i<^' the spin operators of each dot. 

In this configuration, spin blockade is manifest for sufficiently small bias voltage 
such that the state with two electrons in the right dot (the state that supports the 
current) is in resonance with the one electron states in each dot. The current is then 
quenched whenever the electrons in each quantum dot have the same spin polarisation 
and the Pauli exclusion principle avoids inter-dot tunnelling [12]. 

Moreover, considering that any excited state in the right quantum dot has an 
energy above the transport window and that inter-dot tunnelling is spin independent, 
a current can flow only through processes involving the doubly occupied singlet in 
the right quantum dot, |5r) and the inter-dot singlet |5o) = (|T,i) — li,T))/"v/2. 
Thus, the occupation of any inter-dot triplet state, |+) = |t,T), |— ) = U,i) and 
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l^o) = (IT) i) + li) T))/\/2 inhibits the transport to the collector [10], unless, as we will 
disciiss below, the singlet and the triplet subspace mix due to any perturbation, as 
for instance an inhomogeneous magnetic field in the sample, which produces different 
Zeeman splittings within each quantum dot. 



2.1. Master equation and full counting statistics 

The non-equilibrium dynamics of a quantum dot system can be described by means 
of the equation of motion for the reduced density operator p, obtained after tracing 
out the reservoirs in the total density operator R: 

p = trnR=-^tTn[H,R]=Cp, (2) 

where C is the Liouvillian acting on the reduced density operator. In matrix notation, 

it reads pi = MijPj, where pj denotes the density operator in vector notation. 

Accordingly, one can define current super-operators J'± which, when acting on the 
reduced density operator, describe the electron tunnelling from the quantum dots to 
the collector and back, thus yielding positive and negative contributions to the current. 
Then, the current can be formulated as the trace of the current operator [18, 19] 

I = etTs iJp) = e trs [{J+ - J-)p]. (3) 

Note that, in the same way as the super-operator £ can be written as a matrix and p 
as a vector, the trace over the system space can be expressed as multiplication with a 
transposed vector in Liouvillian space, wj, which is the unit matrix in vector notation. 
Thus, the trace condition tr^ p = I reads v^p = 1, while trace conservation tr^ p = 
corresponds to the relation 

vie = (4) 

and the current expectation value becomes / = ev\jp. 

It is convenient to write the master equation with the help of the current 
operators, which is easily done by identifying those terms that change the number 
of particles in the collector included in J±. Then the master equation assumes the 
form 

p = C{t)p={Lo{t) + J+ + J-)p. (5) 

In a quantum dot system, the super-operator £o(t) describes both the electron 
dynamics and tunnelling through the emitter barrier. 

Both current and shot noise can be expressed in terms of the accumulated charge 
in the collector, eN{t), such that [20] 

/=e^(7V(t)), (6) 

S = e'^^{{N^{t))-{N{t)Y). (7) 

Accordingly, one can obtain expressions for the higher order moments of the current 
by evaluating the expectation value (A'"") that define the statistics of the transmitted 
electrons — the full counting statistics [21,22]. 

For the specific computation of the cumulants, we define the operator 

G{z,t)=irn{z''R{t)), (8) 



Spin correlations in spin blockade 



4 



which generahses the density operator. The latter is recovered in the hmit G{z —> 
l,t) = tr'jiR{t) = p{t). G{z,t) obeys the equation of motion 

G(z, t) = (£ + (z - 1)J+ + (z-i - 1) J^) G{z, t), (9) 

while its trace is the moment generating function for the transported electrons, which 



where 



g(")W= (z^) G{z,t) 



trsg^'^Ht), (10) 



tr^ (iV"i?(i)) . (11) 



Introducing the shorthand notation g = g^^"^ and g' = g^^\ we find the equations of 
motion [19] 

m = p{t) = cp{t) (12) 

g'{t) = Cg\t) + {J+-J_)p{r). (13) 

The corresponding equations of motion for the higher-order moments allow the 
recursive computation of the full counting statistics without the need of explicitly 
computing eigenfunctions and derivatives [23]. This is particularly helpful for larger 
systems that require a numerical treatment. 

2.2. Stationary solution and zero frequency noise 

The average current and the shot noise can be obtained by integrating the equations 
of motion (fT2| and (fT3|) . For long times, the system evolves into a stationary state 
determined by Cpoo = 0. The stationary solution poa thus is the eigenvector of the 
Liouvillian that corresponds to the eigenvalue zero with the normalisation v^poo = 1- 
By inserting this into equation p3p . we obtain 

g'{t)^Cg\t) + {J+-J-)poo. (14) 

The eigenvalue zero of £, which corresponds to the stationary state, involves a solution 
with one component linear in time, which can be singled out by projection to the null- 
space, PooVq' such that 

g'{t) = Poovl{J+ - J-)poot + g'At)- (15) 

The orthogonal component g'^t), converges for long times, as we will demonstrate 
below. The information about the shot noise is fully contained in g'j_{t), and for its 
stationary solution, we obtain by inserting (|15p into (jl3p the algebraic equation 

Cg'Acx^) = ipoovl-l)iJ+^J-)poo, (16) 

together with the orthogonality condition uJ(7^(oo) = 0. The numerical integration 
of equations (|12p and as well as the projection to the orthogonal subspace are 
illustrated in figure [H 

Knowing paa and g'^oo) allows one to compute the stationary current and the 
zero-frequency noise 

I = e vl{J+ - J^Poo (17) 
S = e^vl{{J+ - J-)g'A^) + {J+ + J-)Poo}, (18) 
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Figure 2. Time evolution of a typical diagonal element of the density matrix p 
and the corresponding terms of g' and g'j_ , see section 12.21 



respectively. The ratio between them defines the Fano factor 
e\I\ 

which reflects the sub- or super-Poissonian character of the noise according to < 1 
or F > 1, respectively. This method is also valid for ac driven systems provided that 
one consideres the averages / and S over one period of the field [19]. 

2. 3. Spin current fluctuations 

Up to now we considered the statistics of the total transferred charge. However, we 
are also interested in the spin degree of freedom an the eventual spin-spin correlations. 
It is straightforward to adapt the formalism presented above to the case of transport 
of electrons with spin a = 1,1 by replacing in equations ([6]) and ([7]) the electron 
number operator, N, by number operator for electrons with spin a = N^- 
Then one obtains the spin current operators J7cr± and the corresponding spin currents 
la = e-^{Na), as weU as the spin auto-correlation Sa = e'^-^{{N^) — {N^)^) with the 
spin Fano factor = 5^/(61/^1). 

Moreover, by comparing charge noise and spin noise, one can also study cross- 
correlations of currents with opposite spin whose zero-frequency component reads 

Sn = ^((^T^i) - (^T>(^^i» = ^(^- E ^-)- (20) 

This quantity allows one to define the dimensionless spin correlation coefficient 

(21) 



which assumes the value ?' = 1 if the transport of "f -electrons and |-electrons is perfectly 
correlated, while r — for uncorrelated spin currents. 

3. Spin current noise for undriven quantum dots 

In the configuration described in section [21 i.e. for two quantum dots with up to two 
electrons in the system, it is convenient to decompose the density operator into the 
basis |1) = |0,T), |2) = |0,i), \a) = |T,i), \b) = |i,T), |5r) = |0,n), |+) = |T,T), 
K) lii i)- As discussed above, transport suffers from spin blockade unless the inter- 
dot triplet states have finite lifetime. However, as discussed in [24], transport becomes 
possible for a small but finite thermal smearing of the Fermi surface of the emitter. 



Spin correlations in spin blockade 



6 



Then an electron in the left quantum dot being in one of the states \a), \b), |±) can 
tunnel with a rate xFl to the emitter, where x ~ 1 — /(cL + t^LR); a-nd can be replaced 
by an electron with opposite spin. We neglect here the contribution of leakage currents 
due to co-tunnelling [25] . 

Under the usual Born-Markov approximation [26], the master equation for the 
reduced density matrix, p = £p, in presence of an homogeneous magnetic field, i.e., 
for Al = Ar reads 



for the off-diagonal elements, i.e. the coherences, where £ = £l — er + C/lr — C^r- 

The spin current operators jT^i can be written as matrices for which all elements 
but (J7|+)2,Sr = {'Ji+)i.Sn = Tr vanish. Processes which transport an electron from 
the collector to the state [Sr) are energetically forbidden and, thus, Ja- = 0. 

Under the condition e = 0, we find that \a) and \b) are in resonance with [Sr) 
such that for symmetric dot-lead coupling Fl = Tr = F, a leakage current / — xeF/3 
emerges. Then, the Fano factor is determined by the internal dynamics even in the 
limit of small thermal excitation a; <i; 1 and reads 



This super-Poissonian shot noise can be explained in terms of electron bunching: Spin 
blockade occurs once electrons decay to a triplet state. Then, the system will remain 
for a rather long period in the triplet subspace, intermitted only by short lapses of time 
in which an electron tunnels to the emitter and is replaced by an electron with opposite 
spin and, thereby, resolves spin blockade. The resulting electron bunching turns into 
super-Poissonian statistics, as in the case of dynamical channel blockade [24, 27, 28] , 
where two interacting channels block each other. Here the two channels are replaced 
by the two interacting spin degrees of freedom. 

Interestingly enough, this result c;an be; c;onipared to the transport through a 
single quantum dot for which one level (the singlet channel) lies within the transport 
window and three levels (the three inter-dot triplets) are below both Fermi surfaces, 
such that the occupation of the later states induces dynamical channel blockade. Then 
the Fano factor becomes F = 1 + (4/3)Fl/(Fl -t- Fr) and for Fl = Fr assumes the 
value F = 5/3. 



Pi =^RPSn+xTi,{pb + p+)-2{l-x)pi 

P2 = TrPSr -I- xTi^iPa + P-) - 2(1 - X)p2 

Pa = - 2tLR ImpSaa - xT^Pa + (1 - x)rj,p2 

Pb = 2<LR Im psjib - xTi^pb + {1- a;)FLPi 

PSr = 2tLR iMPSna - PS^b) - ^^RPSb. 

p+ = (1 - a;)FLpi - xTlp+ 

P- = (1 - X)T-Lp2 - xTi^p- 



(22) 



for the diagonal elements which describe occupation probabilities and 

Pab = iiLR(pSR6 + PaSn) - xTj^pab 





(24) 
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Figure 3. Current (in mV/?i), Fano factors and spin correlation for the undriven 
case, as a function of Zeeman inliomogeneity for (lr = 0.0026, el = 1.5, 
£pi = 0.45, Uji = 1.45, Ar = 0.026 (corresponding to a magnetic field of ~ IT), 
/iL = 1-94, /.tR = 1.1, Fl = Fr = 10-=* and k-^T = 0.001 (energies in meV). 



In the absence of spin scattering processes, both spins contribute equally to the 
current, = //2, and the Fano factor for each spin current reads 

Moreover, we obtain the spin correlation coefficient 

^x + 0{x^). (26) 



As can be seen in figure [31 this behaviour considerably changes as soon as an 
inhomogeneity in the Zeeman splittings appears. Then mixing between |S'o) and \Tq) 
adds a new conducting channel while the number of blocking states is reduced to the 
two states |+) and |— ) and, thus, the effective transmitted charge and the Fano factor 
increase. Then we find F = 2>, F„ = 2, and r = 1/2. 



4. Electron spin resonance 



The spin of an electron can be manipulated by external magnetic fields. For instance, 
in the presence of an oscillating magnetic field Sac whose frequency equals the Zeeman 
splitting produced by a constant magnetic field, the electron spins rotate, which is 
know as electron spin resonance [29]. For quantum dots, this effect may manifest itself 
in current oscillations [8] and allows one to determine the spin scattering times [30]. 

In a double quantum dot, one expects that the rotation of a single electron 
spin removes spin blockade [9], but the presence of a second electron can lead to 
collective rotations which quench the current or affect the current oscillations [10]. 
This behaviour is reflected also in the shot noise characteristics, as we will find below. 

When a magnetic field with an ac component circularly polarised perpendicular 
to the x-y-plane acts, then the z-components of the electron spins start to rotate, as 
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Figure 4. Spin rotation in a double quantum dot induced by an oscillating 
magnetic field at resonance with the Zceman splitting. 



depicted figure ID The corresponding time-dependent Hamiltonian reads 

i?B [t) = ^ [a, 5: + ffSac (Sl COS cut + SI sin c^i) ] . (27) 

i 

In order to rotate the electron spin, the frequency of the ac magnetic field must 
satisfy the resonance condition huj = A; [8]. Then one would naively expect that 
the Pauli exclusion principle eventually would not apply to the interdot tunnelling 
such that spatial oscillations between the left and right quantum dot occurred and 
a finite current emerged. For a homogeneous Zeeman splitting through the double 
dot, Al = Ar, however, this is not the case. As the oscillating magnetic field 
brings different spin projections within one dot into resonance, the electron spins 
in different quantum dots rotate simultaneously within the triplet subspace, which is 
decoupled from the singlets. Then, the electrons become eventually trapped in the 
triplet subspace and the tunnelling current towards the collector drops to zero [10]. 

The Zeeman splittings in different dots may as well be different, Al ^ Ar, for 
example due to an inhomogeneous magnetic field B^c- A further reason may be a 
dependence of the g- factors on the particular quantum dot, or a difference in the 
hyperfine interaction [13,31]. Then, the states \So) and \Tq) mix and, thus, only 
the triplet states 1+) and |— ) suffer from spin blockade. Moreover, since for an 
inhomogeneous Zeeman splitting, _Bac can only be at resonance with the electrons 
in one of the quantum dots, the trapping in the triplet subspace is lifted and a finite 
current can flow; see figure [S] 

The contribution of these two effects — single electron spin resonance and the 
suspension of channel blocking by singlet-triplet mixing — results in a non-monotonic 
dependence of the Fano factor, which is considerably enhanced in the vicinity of the 
degenerate Zeeman splitting Al ~ Ar . In that region, the rotation of the two electrons 
still dominates and increases the blocking time which is much longer than the time 
lapses with finite conduction. This leads to an enhanced electron bunching which is 
manifest in the larger Fano factor shown in figure [5] This super-Poissonian behaviour 
is accompanied by an enhanced spin correlation coefficient r. Once single-electron 
spin resonance dominates owing to large Zeeman inhomogeneities, the Fano factors of 
both the charge current and the spin current saturate at the values F = 3 and = 2, 
respectively, while the spin correlation coefficient becomes r = 1/2. 



4.I. Crossover to the incoherent regime: unblocking by relaxation 

The presence of inelastic spin scattering processes damps the coherent rotations 
within the triplet subspace resulting in a finite population of the singlet subspace. 
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Figure 5. Current (in mV/?i), Fano factors, and spin correlation as a function 
of the Zccman inhomogeneity for magnetic field intensity Bac = 0.1\/2T, iLjs = 
0.0026, El = 1-5, £r = 0.45, C/r = 1.45, hu) = Ar = 0.026, /.jl = 2, ^iR = 1.1, 
Tl = Tr = 10-3, g = 0.026 (energies in meV). 



Then, a finite current emerges. We treat relaxation processes phenomenologically by 
introducing in the master equation spin relaxation and spin dephasing times, Ti and 
T2, respectively. 

The results for the Fano factors and the spin correlation for Al — Ar shown 
in figure [6] indicate the existence of two different regimes: When spin relaxation 
processes represent merely a perturbation, only a small leakage current flows and the 
Fano factors approach = 5/3 and F^, — 4/3, while the spin correlation coefficient 
becomes r = 1/4. As spin relaxation becomes more important, the current increases to 
a maximum which marks the crossover to the incoherent regime. At the maximum, the 
current and the spin current are sub-Poissonian, F^F„ < 1. For even faster relaxation 
processes, typically for Ti < 27r/r2ac, both spin rotation and inter-dot tunnelling — 
which are considered to be of the same order here, with Rabi frequencies fiac — 2i3ac 
and JIlr = 2V2^lr, [10] — are no longer effective, so that spin relaxation takes place 
immediately after spin rotation. Consequently, no electron flow is observed [10]. As 
expected, in this regime the currents are almost Poissonian and the spin projections 
are uncorrelated. 

5. Pumping and photon-assisted transport 

Not only an ac magnetic fleld, but also an ac electric field can support the electron 
transport by exciting electrons to states which lie above the chemical potentials of 
both leads and, thus, are energetically forbidden. Then, the triplet states have finite 
lifetimes due to photon absorption processes so that spin blockade is resolved [1,32]. 
A paradigmatic setup for which resonant ac electric fields can lead to finite charge [33] 
or spin currents [16] even in the absence of a source-drain voltage (/^l = = m) is 
the spatially asymmetric pump configuration of the double dots sketched in figure [T] 
There the asymmetry is provided by the different intra-dot interactions C/l > U^. 
Pumping occurs when the frequency of the electric ac fleld provides the energy 
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Figure 6. Current (in mV//i), Fano factors and spin correlation as a function of 
the spin relaxation rate. We have considered T2 = O.lTi and the same parameters 
as in figure |5] except for Al = Aj^ = 0.026. 



/"l Ilr Tr 




Figure 7. Pumping configuration, where spin blockade is removed by photon 
assisted processes through the contact barriers. 



necessary to transfer the electron from the left quantum dot to an energetically higher 
singlet state |S'r) in the right dot. For this process, the resonance condition reads 
huj = TiLOfi = — £l + Utx, — ?7lr- As before, the singly occupied states of both dots 
are below both chemical potentials, e; < /.t, while the doubly occupied singlets lie well 
above, [// + £; > /i. We restrict ourselves to the pumping configuration, because 
increasing the bias voltage reduces the contribution of photon assisted processes 
through the contact barriers [34,35] which are essential for the removing spin blockade. 

The ac electric field is modelled as an oscillatory dipole potential, i.e. as a time- 
dependent energy shift with a phase tt between the two dots. Then the Hamiltonian 
([1]) acquires a term 

^ac(0 = ^("L-nR)cOs(wt). (28) 

By the unitary transformation U{t) = e'('=^'"^/2?ic^)(riL-nR)sini^t^ ^Yi^ ^j^^g dependence 
is transferred to the tunnel couplings according to 

H'{t) = U{t) (H^lhdt) UHt) = Ho + H'LR{t)+H'Tit) + Hle.dsA'2^) 

where 

i?'LRW= E (-l)'''/J^)E(*LR^'''"*^LcR. + H.c.), (30) 
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Thus the driving effectively renormalises the inter-dot tunnelling by the Bessel function 
J^, where the index v reflects the number of photons involved. This creates dynamical 
charge localisation under conditions that we specify below. The corresponding photon- 
assisted tunnelling rates turn out to be 



where the former rate governs processes \n) — > |m) that remove an electron from the 
double dot, while the latter refers to adding an electron to the system. The density 
of states in the leads, di, is assumed to be constant. 

If a magnetic field is applied such that the energies of the spin-down electrons 
are shifted by A;, the spin-up electron is delocalised within the double quantum dot 
occupying the states |t,i) and |0, ti) until one of the electrons tunnels to the right 
lead. If the Zeeman splitting is the same in both dots, the spin-down electron can at 
the same time be delocalised between |0, Sr) and ||, t) for the same frequency. 

At low ac intensities, photon-assisted tunnelling is not significant. Indeed, using 
the approximation Jv(x) ~ 5^,0 for small x, we find that photon-assisted tunnelling 
rates become identical to the rates for the static system discussed in section [S] cf. 
figure m When increasing the intensity, both the current and the noise exhibit a more 
involved behaviour stemming from the amplitude and frequency dependence of the 
tunnel rates via the arguments of the Bessel functions Ji{eVac/?icu) and Ji(eI4c/2?iw) 
which govern inter-dot tunnelling and photon-assisted tunnelling, respectively. When 
the ac intensity approaches the value for which Ji(eVac/?iw) assumes its maximum, 
inter-dot tunnelling is most effective and, thus, the current assumes its maximum 
as well, while the Fano factor assumes a minimum [36]. The spin Fano factor, 
by contrast, does not assume a minimum and, consequently, the spin correlation r 
becomes negative. Increasing the ac intensity further reduces the net current because 
in this limit, the rates of left-to-right tunnelling approach those of the opposite 
processes. This cancellation does not affect the noise strength S, so that a smaller 
current leads to an increasing Fano factor. 

Particularly interesting is the condition Ji(eVac/fu^) — for which dynamical 
charge localisation is found. Then one-photon interdot tunnelling is suppressed due 
to the renormalisation factor in the hopping term (j30p . Since for our configuration 
these processes represent the main contribution to the transport, the current becomes 
strongly suppressed, as can be seen in figure [5] in the region close to eT4c ~ 3.8huj. 
Nevertheless, photon-assisted tunnelling through the contacts is still noticeable by its 
contribution to the current fluctuations, leading to a sharp peak in the Fano factor [37]. 
This behaviour is different from the one found in the large bias limit, where current 
suppression is associated with Poissonian noise [38] . 

Another important factor for spin blockade is the difference between the Zeeman 
splittings, as discussed in section H) If Al = Ar, the states |t,i) and ||,t) are 
indistinguishable, then the inter-dot singlet \So) and the triplet |To) influence the 
dynamics. Since inter-dot tunnelling does not change the total spin, it can only occur 



r 



r. 



mn 



mn 




(1 - Mhiv^n + lyhw)) , (32) 

\jl\'^fl{hLL>mn + l'h(^), (33) 
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Figure 8. Current (in mV/?i), Fano factors and spin correlation as a function of 
the electric field intensity for u) = ljq, ^lr = 0.005, £l = 0.4, ej^ = 0.25, Un = 1, 
(7lh = 0.5, Al = Ar = 0.026, ml = Mr = 1-2, Tl = Tr = 10-=^ (energies in 
meV). 



within the singlet states l^o) and |S'r). Thus, not only the states |t,T) and |J,,|) 
contribute to the transport blocking, but also the triplet state ITq)- 

On the other hand, if the Zeeman splittings are different, Al ^ Ar, spin blockade 
is less effective due to the mixing between the states ITq) and \Sq) [10]. Then the 
current increases, as can be seen in figure [51 The reduced the number of blocking 
states again diminishes the Fano factor. If the difference between the Zeeman splittings 
is large enough, we find resonant one-photon inter-dot tunnelling for electrons with 
one particular spin polarisation — here the spin-up polarisation. Consequently, the 
transport becomes spin dependent, as can be appreciated in the different Fano factors 
in figure [5) 

6. Conclusions 

In this work we have studied the suspension of spin blockade in the electron transport 
through coherently coupled double quantum dots by ac driving fields and also 
by thermal excitations. In particular, we focussed on the associated shot noise 
characterised by the Fano factors for the electron and the spin currents, as well as 
on the spin-spin correlation coefficient. 

One way to resolve spin blockade is the application of magnetic fields that cause 
Zeeman splitting and spin rotation. Most interestingly, we found that the current 
noise for both charge transport and spin transport depends sensitively on whether the 
Zeeman splitting is the same or different for the two dots. For degenerate Zeeman 
splitting, the eigenstates are delocalised and the transport across the whole sample is 
dominated by tunnel events which obey Poissonian statistics. Thus the Fano factors 
tend to be close to unity. If the degeneracy is absent, transport is dominated by lapses 
of time at which spin blockade is lifted. During these lapses, we observe the transport 
of bunches of electrons with correlated spins. 

Considering also spin relaxation, we find that such incoherent processes can even 
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Figure 9. Current, Fano factor and spin correlation as functions of tlie Zeeman 
inhomogeneity (by varying Al) for eVkc = fti^ = Ty-^o ^.nd the same parameters as 
in figure [8] Though the spin currents almost coincide, as seen in the upper panel, 
their dynamics can be distinguished in the noise components. 



contribute to the suspension of spin blockade: With an increasing relaxation rate, 
transport is enhanced and becomes more regular and even sub-Poissonian. When the 
predominantly incoherent regime is entered, the Fano factor increases again until we 
eventually reach an incoherent regime with Poissonian tunnel currents. 

Time-dependent fields can also resonantly excite electrons to orbitals with energies 
above both Fermi surfaces. In that way, they may cause photon-assisted transport and, 
in asymmetric configurations, electron pumping. The latter results in a contribution 
to the dc current at zero bias voltage. Photon assisted transport can be described by 
performing a unitary transformation of the Hamiltonian such that only the tunnelling 
matrix elements are explicitly time dependent. This renormalises the tunnel couplings 
such that they depend on driving amplitude and frequency. The main contribution to 
interdot hopping can even vanish, so that for specific driving parameters, the electrons 
become localised within one dot. Consequently, the current is almost suppressed. 
Nevertheless, the current fluctuations may stay at a significant level and, accordingly, 
the shot noise level is super-Poissonian. With an additional Zeeman splitting, the 
resonance condition becomes spin dependent. Then photon- assisted transport favours 
one particular spin projection and the two spin currents are no longer identical. This 
also affects the spin shot noise, although the difference in the corresponding Fano 
factors turns out to be relatively small. 

In conclusion, our results underline that ac fields can have intriguing consequences 
for spin-dependent transport. Depending on the type of excitation, both charge and 
spin transport can occur in bunches. We have shown that this should be clearly visible 
in the corresponding measurements of the Fano factor and, most likely, also in higher 
orders of the full-counting statistics. For driven spin transport this still is awaiting 
closer investigation. 
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